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ABSTRACT 


A solution is presented to the radiative transfer of the .?olar 
irradiation through a turbid atmosphere, based on the single- 
scattering approximation, i. e. , an assumption that a photon 
that underwent scattering either leaves the top of the atmos- 
phere or strikes the surface. The solution depends on a 
special idealization of the scattering phase function of the 
aerosols. The equations developed here are subsequently 
applied to analyze quantitatively (1) the enhancement of the 
surface irradiation and (2) the enhancement of the scattered 
radiant emittance as seen from above the atmosphere, caused 
by the surface reflectance and atmospheric back-scattering. 
An order-of-magnitude error analysis is presented. 
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FIGURE CAPTIONS 


Figure 1. The increase in the surface irradiance by surface reflection and 
subsequent backscattering, Sj.||j, for (i) Q = 0. 1, f/b = 1. 0; (ii) Q = 0.3, 
f/b - 3. 0; (iii) Q = 0.3, f/b = 5. 0. 0 = 0° 

Figure 2. The increase in the surface irradiance by surface reflection and 
subsequent backscattering, for (i) Q = 0. 1, f/b = 1. 0; (ii) Q = 0. 3, 
f/b = 3. 0; (iii) Q = 0. 3, f/b = 5. 0. 0 = 60° 
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Figure 3. The increase in the atmospheric radiant emittance by surface 
reflection and subsequent upward scattering, Srf* for (1) Q * 0. 1, f/b » 1. 0 
(ii) Q = 0. 3, f/b = 3. 0; (iii) Q = 0. 3, f/b = 5. 0. 0=0® 

Figure 4. The increase in the atmospheric radiant emittance by surface 
reflection and subsequent upward scattering, Srf, for (i) Q= 0.1, f/b= 1.0 
(ii) Q = 0. 3, f/b = 3. 0; (iii) Q = 0. 3, f/b = 5. 0. 0 = 60® 
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SINGLE SCATTERING SOLUTION FOR RADIATIVE TRANSFER 
THROUGH RAYLEIGH AND AEROSOL ATMOSPHERE 


INTRODUCTION 

The problem of radiative transfer through the Earth* s atmosphere is highly 
complex - and an exact solution defies a simple analytical treatment. For a 
survey of approximation methods see Hunt (1971). In what follows we present 
a simple approximate solution for solar irradiation based on single- scattering, 
i. e. , an assumption that a photon once it is scattered, either leaves the at- 
mosphere or strikes the surface. At the surface it can be either absorbed 
or reflected. The solution can be regarded as the limiting case for an optically 
thin atmosphere. Furthermore, the simplicity of the solution hinges on an 
idealization of the aerosol scattering phase function, which is first described. 

The equations derived here are applied to analyze the increase in the surface 
irradiance due to a reflection from a bright surface and subsequent back- 
scattering, and also to quantitatively describe the inhomogeneities in the atmos- 
pheric veil introduced by surface reflection and subsequent upward scattering. 
These inhomogeneities will be analyzed more fully in a forthcoming paper 
(Otterman and Fraser, 1978). Also detailed treatmert of the atmospheric ef- 
fects on imaging in Landsat bands is planned (Fraser, Kaufmant Otterman and 
Podolak, 1978). 

SCATTERING PHASE FUNCTIONS FOR RAYLEIGH AND 
AEROSOL SCATTERING 

The angular dependence of scattering of unpolarized light by the air molecules 
is described by the well known phase function 1 + cos^ 0 , where the angle 0 
is measured from the incident beam. Integration of this expression over a 

sphere, i. e. , forming the integral 2 ir f sin^ (1 + cos^ 0 ) d 0 = 47 t + ~ . 

o 3 

indicates easily that in Rayleigh scattering % of the scattering is isotropic and 
only V 4 predominantly forward or backward. 

The phase function of aerosol scattering can vary greatly depending on their 
electro- magnetic properties, on the size distribution of the particles, and 
their shape (spherical or more needle like) (Kondratyev, 1969). The phase 
function is here greatly idealized, and is represented by four parameters a , 

|3 , 7 , m as follows: 
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i The fraction a of the scattered photons is limited to a very 
narrow forward angle. 

ii The fraction /? is scattered in direction 0 , into a eolid angle 
2 JT sin 0 d 0, proportionally to the product of this angle and of 
cos^0 (m is even). 

iii The fraction 7 is scattered isotropically. 

Since a + + 7 = 1, only three of the parameters (including m) are in- 

dependent, and have to be determined appropriately by fitting to a particular 
phase function. 

In this idealized phase function the ratio of forward scattering to back-scattc. *- 
ing by the aerosols (i. e. , the ratio of scattering into the hemisphere O<0<?r/2 
to scattering into the hemisphere Tr/2<0<ir; in our model also the ratio of 
scattering into any hemisphere that includes 0 = 0 to the scattering into the 
opposite hemisohere) is given by (1 + a )/(i - a). 


DERIVATION OF THE ANALYTICAL EXPRESSIONS 


The solar irradiance on a horizontal surface at the top of the atmosphere is 
given by cos 0 ^, where is the solar zenith angle. A plane-parallel atmos- 
phere is assumed; thus the airmass is given by l/cos 6 ^. Single scattering is 
assumed, that is, a photon that is scattered, without undergoing a second en- 
counter impinges on the ground or leaves the atmosphere. The direct irradiance 
of the surface, G^, is 

-Q/cos 

G^j = cos 0Q e (1) 


where Q is the vertical optical thickness at the wavelength considered. 

The components of the optical thickness are Rayleigh scattering, R, aerosol 
scattering, M, and absorption B: 


Q = R + M + B (2) 

According to our idealized phase function, the fraction [ R + (l + a)M ] /2Q 
of the photons that are scattered or absorbed from the direct beam, are scat- 
tered into a hemisphere below a horizontal plane through the scattering point. 
Because of the single scattering assumption, this fraction reaches the surface. 
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Therefore the irradiance Gg^ contributed by the scattering from the direct 
beam to the surface is: 


Ggd= cos^o 


/ -Q/cos 6 oV 

r ) 


2Q 


( 3 ) 


The Irradiance of the surface, G^, when no surface reflection is considered, 
is: 


G = Gj + 
^o d 


°sd= 


-Q/cO 50 ^j 

e 




a+g)U 

2Q 


(4) 


A fraction of the irradiance of the surface is reflected, and this fraction in 
part is back- scattered to the surface. This reflection and back-scattering 
from the atmosphere form a geometric series effect: 


+ kG„ + k2G„ ^ 


G 


o 


1 - k 


(5) 


where G^ is the total irradiance and k the coefficient of the reflection back- 
scattering. This coefficient, assuming reflectance according to the Lambert's 
law, with reflectivity a^, is given by: 


it/2 / -Q/COS0 \ 

_ ^ J COS 4> 8in0 - e J 


(l-a)M 

2Q 


dtp 


R+ (l-a)M ^ 

g *-'1 

O Q 


( 6 ) 


since [ R + (1 - a ) M 1 /2Q is the fraction of the scattered or absorbed photons 
that is back-scattered downward, and where 


/ -Q/cosa\ 

(Q) = f cos^ sin^ll-e Jd0 

^ ' ** n n 

= ‘/z 1 - (1 - Q) + q 2 (+0.5772 + InQ + ) 

L \ n=i /J 


(7) 
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The function Cj(Q) is tabulated In Appendix A. For Q < i, Cj(Q) ^ Q. 
Thus the total irradiance is:* 


G. = cos 
t o 


-Q/cos(?o / -Q/COS0 \ (1+ g) 

+ \l-e / 2Q 

l-^o C,(Q) 


M 


(5. 1) 


which for a Rayleigh atmosphere, Q = R, can be simplified to 


Gt= COS0Q 


-R/cos0. 

1+ e ° 


2[l-a^Ci(R)l 


(5. ?) 


The absorption Eq at the surface is given by 


Eo- (1-a^) G^ 
and in the atmosphere, E„ 

a 

7t/2 


( 8 ) 


2 7ra„G, 




^ j cos« sin* ^ _ ^-Q/cos9„ j ^ 


2 GtaoB 

Q 


Ci(Q) + cos 




( 9 ) 


We analyze the regional Earth-atmosphere system hemispheric reflectivity 
A^ in terms of the radiant emittance from the top of the atmosphere A^cosAq, 
which is given by the irradiance cos less the absorption terms: 


Atcoafl^ = cosOq - (Eq + E^) 


cos 0, 


1- (1-e 


•Q/cos0,A ^ 

) Q 


2aoB 


1-ao *■ 


Q 


Ci(Q) 


-Q/cosO 




-Q/cosO \ R+ (l + g)M 




1 - a R^ (1-Qfi . M Cj(Q) 
° Q 


( 10 ) 
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The planetary albedo, A-, for the illuminated hemisphere (i. e, 0^=0 is 
located at the center of me hemisphere) is: 


ff /2 


2ir Jsin0QAjCO8 5^d0Q 




■wjr 




sln0«cos0^ dfl, 


2a^B 


2 _ ji + — 

o Q 


C,(Q) 


1 - 2C,(Q) + C,(Q) 


S. 




2B 

Q Cj(Q) (11) 


where the term (1 -a^ 

2a B 
o_ 

Q 


Cj(Q) X . . . + 


) X. 

2B 

Q 


. . . represents the absorption at the surface and 
Ci(Q) represents the absorption in the atmosphere. 


The radiance at nadir (for an arbitrary consists of three elements; the 
"signal" radiance Ljj j. from the "target" pixel with reflectivity r; 


rG^e 


n r 


( 12 ) 


the radiance Ljj scattered from the direct beam 


■-Q/COS0J 


^nd= cos0o(l-e 


((m+l)^cos™0^+7]^+ (cos20^^1)y||- 


(i;j) 


and the radiance ^ reflected from the terrain In the general vicinity of the 
target (reflectivity a) , and subsequently scattered to nadir; 


Gja 


"n a 


X 


^l-e"^ja^ + 27 t I j sin0 cos*^ 

o 

{p(m*D cos”>« * t)^ * (cos“« * I) 


X 
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( 14 ) 



P{m+ l)Cm+i (Q) + 7C,(Q) 



+ 


Cg(Q) + C,(Q) 


8Q 


where 


nl2 

Cn,+i(Q) = I d0 

0 

These functions C^+j are discussed in Appendix A. 

The radiance at nadir is 


(15) 


L 

n 


= L. 


n r 


+ L 


h d 


+ L 


n a 


(16) 


APPUCATIONS 


Enhancement of Irradiance of the Horizontal Surface by the Surface 
Reflectance and Subsequent Back-Scattering 

In what follows we introduce two abbreviations: the fraction f of the atmos- 
pheric optical thickness that effects forward scattering 

2Q ' 


and the fraction b that effects back- scattering 


R^^ a-g)M 
2Q 


(18) 


and thus f/b denotes forward to back- scattering ratio of the atmosphere, and 
1 - (f + b) the absorbing fraction. 
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The surface Irradiance due to scattering from the direct beam is 


G , = cos 6 
sd 




-Q/COJ0Q \ R+ g + g)M _ ^ -Q/cos0^ 


2Q 


iOsSjj ( 1-e 


/ -Q/cose^v 

(... ) 


( 19 ) 


The surface irradiance due to surface reflectance and back* scattering is 


Orb’ 


-<3/cmS„ 




,b C,(Q) 


- 1 


cosa, - (l-e-^^"^o)f| 2a,b C,(Q) 

l“2aob C,(Q) 


( 20 ) 


The enhancement Sj.^j of the skylight by this effect relative to the scattr ■ 
from the direct beam therefore is 


^rb G 


’£b 

sd 


Sob |e 


-Q/cos0p / -O/COS0, 


f|2C.(Q, _ 


f(l _ e^Q/cos^o^ jl-2aobCi(Q) 


a b 
o 


*Q/cos0„ 


1 + 




2 C<(Q) 


l-2a b C,(Q) 


( 21 ) 


It can be seen that the enhancement is proportional to a^b, but for low Q and 
small zenith angles 0^ (high solar elevations) it is proportional to a^b/f , 

because then — - — — — — ^ IL > 1 . 

Under these conditions (since C](Q) % Q for low Q) 


7 



and is thus approximately independent of Q. 

The fur ction Eq. 21, is plotted vs. for - o“ (see Figure 1) 

i) for Q = 0, 1, and for f = b = 1/2, i, e. , for a Rayleigh atmosphere about 
0.55 fim 

ii) for Q = 0. 3 and for f = 3/4, b = 1/4, which for scattering by aerosols 
only corresponds to oi~ 0. 5; or when Rayleigh optical thickness is 0. 1 
and aerosol optical thickness is 0.2, corresponds to a = 0. 75, 

iii) for Q = 0. 3 for f = 5/6, b = 1/6, which for scattering by aerosols only 
corresponds to a = 2/3; or when Rayleigh optical thickness is 0. 1 and 
aerosol optical thickness is 0.2, corresponds to a = 1.0. 

In Figure 2 the same graphs are repeated for 0® = 60®. 

Coulson (1968) pointed out that many mineral surfaces have high reflectivity 
at the longer visible and solar Infrared wavelengths, and over such surfaces, 
surface reflection and subsequent back-scattering from a Rayleigh atmosphere 
can contribute as much as half of the total skylight particularly at high sun 
elevations. Indeed, in Figure 1 it can be seen that for a Rayleigh atmosphere 
one-half of the total skylight (Sj.jj = 1. 0) is contributed by this effect at a^ = 
0.57. Such reflectivity can be approached in sandy desert areas at wavelengths 
0. 7 - 1. 1 fitn (Otterman and Fraser, 1976), can be encountered over thickly 
vegetated terrain above 0.75 Atm (Colwell, 1974) and can be exceeded over 
snow at any wavelength. For a real atmosphere, i.e. Rayleigh and aerosol, 
a reduction of this effect occurs in the ratio b/f, see Figure 1. For a higher 
zenith angle the effect is less significant relatively to the scattering from the 
direct beam, see Figure 2, but remains significant in terms of contribution to 
the total irradiance, 1. e, , as compared to the direct irradlance. 


Enhancement of the Atmospheric Veil by the Surface Reflectance and 
Subsequent Scattering 

The enhancement of the atmospheric veil by the surface reflectance and subse- 
quent scattering forms a complimentary question to the problem just discussed. 
Earlier, we analyzed the fraction of the reflected and scattered radiation that 
impinges on the ground, here we discuss the fraction that escapes from the 


atmosphere. This part of the atmospheric veil appears most strongly above 
bright terrain, and thus spatial inhom(^eneities in this atmospheric veil are 
introduced even if the atmosphere is horizontally homogeneous. &ich in- 
homogeneities can introduce difficulties in thematic mapping by multispectral 
radiometry approach and this aspect requires more detailed treatment, in 
terms of radiances measured from a satellite. Here we analyze the enhance- 
ment in the radiant emittance emerging from the top of the atmosphere above 
a bright terrain. 

The radiant emittance due to the scattering from the direct beam is 

„ /. -Q/cos0„\ R+(l-a)M 

"sd = °) 

“ cos9^ ^ 1 jb (23) 

and the radiant emittance from surface reflection and subsequent scattering is 


R + (1 + a ) M 


-Q/COS0 

a^cos0 |e 
o o 


/ -O/cos0_\ 1 

(l-e «)f|2fCi(Q) 


1 2aQb Cj(Q) 


(24) 


and the enhancement over bright terrain with reflectivity a^ as compared 
to terrain a^ = o is 


^rf~ VV 


sd 


a f 
o 


-O/COS0 

e + 




(Ci) 


b {l-2a^bCi(Q)j (l-e 




-Q/cos0^ 


1 + 


f(l_e-Q/-^o) 


2 C,(Q) 


l-2a b Ci(Q) 


(2&) 


There is a strong similarity between this expression for S , and that for 
Eq. (21), which is not surprising in view of what was said before. For low Q 
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and small zenith angles (high solar elevations), the enhancement is approxi- 
mately 

2a f 

(26) 

We compute Sj.£, Eq. (25), as a function for a^ for 0® = 0 (see Figure 3) and, 
as before in Figure 1; 

i) for Q = 0. 1, and for f = b = 1/2, i. e. , for a Rayleigh atmosphere 
about 0. 55 pm 

ii) for Q = 0. 3 and for f = 3/4, b = 1/4 

iii) for Q = 0. 3 and for f = 5/6, b = 1/6 

and repeat the calculations for 0°= 60° (see Figure 4). The case (i) Is 
numerically equal in Figure 3 to that in Figure 1, and in Figure 4 to that in 
Figure 2, but case (ii) and (iii) are higher in the ratio f2/b2in the corre- 
sponding figures. This ratio is 9 for (ii) and 25 for (iii). 

It can be seen from Figure 3 that the atmospheric veil can be more than 
doubled over high reflectivity terrain, since Sj-f reaches 1. 0 at values of 
ranging from 0.14 (case ii) to 0.22 (case iii). 

For thematic mapping by multispectral radiometry, the significant question is 
how to expurge the radiance in Eq. (16), which depends on the surface re- 
flectance from the area around a pixel, and to measure the reflectivity r of the 
pixel (or re"^) from Ljj. This will be discussed by Otterman and Fraser (1978). 

DISCUSSION AND CONCLUSIONS 


Analytical expressions have been derived for the radiative transfer through the 
atmosphere for a limiting case of a thin atmosphere. The solution is based on 
a single encounter (scattering or absorption) approximation to the atmospheric 
passage of the solar radiation, both for the solar irradiance from the top of 
the atmosphere downward and for the surface reflection from the bottom of the 
atmosphere upwards. The simplicity of the solution also depends on the ideal- 
ization of the scattering phase function. 

Presentation of such a simplified solution has to be accompanied by an assess- 
ment of its accuracy, and this is outlined here. The "subsequent encounter" 
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(i. e., either an absorption- encounter or a scattering-encounter of a photon 
that first underwent scattering) are analyzed in an order-of- magnitude 
approximation. Consider that all the first encounters as computed here 
occur on the average mid- way through the atmosphere. The once 

/ -Q/cosfl \ 

scattered photons, (a fraction ^1-e M)/Q (R+ M)/cos 6^ of the 

photons impinging on the top of the atmosphere), face^ if they are scattered 
predominantly forward or backward, an absorbing optical thickness of cos 6^, 
The once-scattered photons that undergo subsequent absorption (before possibly 
impinging on the surface) constitute a fraction ~ (R + M) B/2 cos^ 9^ of the 
impinging radiation, as compared to a fraction B/cos that is absorbed 
before undergoing scattering (and which is accurately computed in our solution). 
The relative error in computing absorption is thus of the order of (R4M)/2 cosd 

The second source of error lies in that second scattering results in a change of 
direction. The effective scattering phase function^ becomes somewhat different 
from that assumed in the computations. Repeating the argument that the second 
encounter for the once-scattered photons is approximately proportional to the 
airmass 1/2 cos 0^^, and that the a fraction of the aerosol scattering represents 
only a small change in the direction that can be neglected here, we estimate 
that the fraction of photons that undergoes a change of direction as a result of 
a second scattering is { R+ (1 - a) M ] /2 cos 6^. In the cases analysed of Q = 

0. 3 and 6q = 60° , figures 2 and 4, this fraction would be between 10 and 15 
percent. The net result would be that an appreciable fraction of the scattered 
photons would be directed at a more nearly random direction relative to the first 
scattering point. The effective phase function then becomes spherically less 
asy metric. The important ratio f/b becomes somewhat reduced, but the results 
for a Rayleigh atmosphere would be scarcely affected. 

There is one limitation of applying the formulas that has to be stressed. The 
above error estimates indicate that the errors increase as l/cos 0^. Beside 
this indication, that the solution becomes quite inaccurate at the low solar 
elevations, the assumption that the fraction o of the aerosol scattering, the 
"a beam, ” is scattered forward into a zero-width spherical angle introduces 
a limitation on the lowest solar elevation for which the formulas can be used. 
According to our idealized phase function, the "a beam” scattered from the 
direct beam will always impinge on the surface. In reality, at the very low 
solar elevations, because of the actual finite angular spread of the ” a beam, ” 
part of the beam will be directed upwards. The formulas should not be thus 
used for 0^ > 70°. 

While this error analysis indicates the need for caution in applying the formulas 
derived here, it is claimed that these simple solutions can provide means 
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of quantitatively assessing the atmospheric scattering effects under low and 
even medium turbidity conditions, especially for low absorptivity. 

The equations are applied to analyze the enhancement of the surface Irradlance 
by surface reflectance and subsequent back-scattering. It is concluded that 
for high surface reflectance of 0.57, one-half of the scattered irradlance can 
be contributed by this effect for a Rayleigh atmosphere. For an aerosol and 
Rayleigh atmosphere, the relative contribution of this effect to the scattered 
irradlance is reduced in the ratio f/b, i. e. , ratio of the forward scattering 
optical thickness to the back-scattering optical thickness. 

The scattered radiant emittance as seen from above the atmosphere is also 
enhanced by surface reflectance and subsequent scattering, and for a Rayleigh 
atmosphere the expressions are numerically the same as in the case of the 
enhancement of the surface irradlance. For a turbid atmosphere, i. e. , Ray- 
leigh + aerosol, the effect is increased proportionally to the f/b ratio. Thus, 
over a bright terrain the atmospheric veil is much more pronounced than over 
a dark terrain. These inhomogeneities constitute an important phenomenon, 
affecting the multispectral radiometry from space, and deserve detailed studies. 
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FOOTNOTES 


*The total irradiance on the surface when =» 1 for 0^=0 and for a non- 
absorbing atmosphere is approximately 1 + [R+ (1 - a) M ] /2, and thus is 
higher than at the top of the atmosphere. This Is due to the fact that some 
photons impinge on the surface more than once. Numerically similar results 
were obtained by Deirmendjian and Sekera (1954). 

ti. e. , a function relating the direction probability for photons emerging from 
the top of the atmosphere or striking the surface, to the direction of the 
direct beam. 


13 


REFERENCES 


Abramowitz, M. and I. A. Stegan, Eds. (1964), Handbook of Mathematical Func- 
tions, U. S. Department of Commerce, National Bureau of Standards Applied 
Mathematics Series 55. 

Colwell, J. E. (1974), Vegetation Cancq^y Reflectance, Remote Sensing of Environ- 
ment, 3, 175-183. 

Coulson, K. L. (1968), Effect of Surface Reflection on the Angular and Spectral 
Distribution of Slgrllght, Journal of the Atmospheric Sciences, 25, 759-770. 

Deirmendjian, D. and L. Sekera (1954), Global Radiation Resulting from Multiple 
Scattering in a Rayleigh Atmosphere, Tellus, VI, 382-398. 

Fraser, R. S., J. Y. Kaufman, J. Otterman and M. Podolak (1978), Atmospheric 
Effects on Imaging in the lANDSAT Bands, in preparation, for submission to 
Applied Optics. 

Hunt, G. E. (1971), A Review of Computational Techniques for Analysing the 
Transfer of Radiation Emerging from a Model Cloudy Atmosphere, J. Quant. 
Spectrosc, Radiat. Transfer, 11, 655-690. 

Kondratyev, K. Ya. (1969), Radiation in the Atmosphere, Academic Press, 

New York. 

Otterman, M. and R. S. Fraser (1976), Earth Atmosphere and Surface Reflectiv- 
ities in Arid Regions from LANDSAT MSS Data, Remote Sensing of Environ- 
ment, 5, 247-266. 

Otterman, J. jnd R. S. Fraser (1978), Adjacency Effects on Imaging by 

Surface Reflection and Atmospheric Scattering, in preparation, for 
submission to Applied Optics. 


14 



1.8 

1.6 

1.4 

1.2 


Figure 1. 
for (i) Q = 








T 



Figure 3. The increase in the atmospheric radiant omittance by surface reflection and subsequent upward 
scattering, for (i) = 0.1, fA) ~ 1.0; (ii) Q = 0.3, f/b - 3.0; (iii) Q = 0.3, fA) = 5.0. 0=0“ 
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Figure 4. The increase in the atmospheric radiant emittance by surface reflection and subsequent upward 
scattering, S^.^, for (i) Q = 0.1, i/b - 1.0; (ii) Q = 0.3, f/b = 3.0; (iii) Q = 0.3, fA> = 5.0. Q = 60“ 



APPENDIX A 


The functions C^^iQ) are integrable for 0<Q<1, m>0 as follows: 

-0/cos^ 

Cj^(Q) = 5 sin ^cos*”^ (1 - e )d^ 

0 

ff/2 -Q/cos^ 

= 5 ^ cos’” ^ sin ^ cos"* d 

0 0 


«» -Qx 
e 


m 


* 1 lx-"* 


dx 


(Al) 


where in the second integral a substitution was made: cos^ =~ and therefore 
dx 

sin\|/dij/ = ’ The second integral is an exponential integral E of order 

m + 2: 


oo -Qx 
e 


^ “STz Em+2«) 


(A2) 


1 X 


which can be expressed in terms of the exponential integral of order 1 by the 
following recursion formula: 


+ 1 (Q) = ~ I (Q) ] . 

m + l'^' m m' ' 


(A3) 


In the particular case for computing Ci(Q): 


E3(Q) = '/2 (1 -Q) e‘^ + — Ei(Q) 


(A4) 


E,(Q) can be represented in terms of a power series in Q: 


n^n 


E,(Q) = -0. 5772 - InQ - L 


(A5) 


A-1 



(Eq. 5. 1. 10, p. 229, in Abramowltz and Stegan, 1964), and thus: 


Ci(Q) = >/2 [1 - (l-Q)e"^ - Q^E,(Q)] . 


(A6) 


= >/2 


1-(1-Q)e“^-Q2 


/ ,^(-l)V\ 

I - 0. 5772 - InQ - L 1 

\ 1 nnl / 


Brief tabulation of this function is given below. 


Table A-1 


R 

Ci(Q) 

0. 00 

. 00000 

0. 05 

. 04508 

0. 10 

.08371 

0. 15 

. 11772 

0.20 

. 14805 

0.25 

. 17532 

0. 30 

. 19996 

0. 35 

.22233 

0.40 

.24271 

0.45 

.26134 

0.50 

.27839 

0.55 

.29405 

0.60 

.30845 

0. 65 

.32171 


A-2 



